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Abstract 



The static cylindrically symmetric solutions of the gravitating Abelian Higgs model form a two 
parameter family. In this paper we give a complete classification of the string-like solutions of this 
system. We show that the parameter plane is composed of two different regions with the following 
characteristics: One region contains the standard asymptotically conic cosmic string solutions 
^-j- , together with a second kind of solutions with Melvin-like asymptotic behavior. The other region 

contains two types of solutions with bounded radial extension. The border between the two regions 
Q\ ' is the curve of maximal angular deficit of 2n. 

Os ' 

O PACS: 11.21. +d, 04.20. Jb, 0440.Nr 

a -1 

1 Introduction 

Cosmic strings jj], || were introduced into cosmology by Kibble Zel'dovich Q and Vilenkin || as 
(linear) topological defects, which may have been formed during phase transitions in the early universe. 
Cosmic strings are considered as possible sources for density perturbations and hence for structure 
formation in the universe. 

The simplest and most common field-theoretical model, which is used in order to describe the 
generation of cosmic strings during a phase transition, is the Abelian Higgs model. This model is known 
to have (magnetic) flux tube solutions || , whose gravitational field is represented by an asymptotically 
conic geometry. These are the local cosmic string solutions. 

This conic space-time is a special case of the general static and cylindrically-symmetric vacuum 
solution of Einstein equations , the so-called Kasner solution: 

ds 2 = (kr) 2a dt 2 ~ (kr) 2c dz 2 - dr 2 - (3 2 {kr) 2{b ^ r 2 d(/> 2 (1.1) 

where k sets the length scale while (3 represents the asymptotic structure, as will be discussed below. 
Notice also that (a, 6, c) must satisfy the Kasner conditions: 

a + b + c = a 2 + b 2 +c 2 = l (1.2) 
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The standard conic cosmic string solution |^|, ^) is characterized by an asymptotic behavior given by 



(1.1) with a = c = 0, 6 = 1 which is evidently locally flat. In this case, the parameter represents a 
conic angular deficit |], ||3], which is also related to the mass distribution of the source. 

A first approximation to the relation between the angular deficit 8<j) = 2-7r(l — (3) and the "inertial 
mass" (per unit length) to, of a local string was found to be ||, O, [l^, p^ |: 

Sep = 8ttGto (1.3) 

Further corrections to ([O]) were calculated in the following papers |lj, [l^, [l(| |T^, . 

In order to fully analyse the string-like solutions of the Abelian Higgs system one needs to solve 
the full system of coupled field equations for the gravitational field and matter fields (scalar + vector) . 
But this is not necessary in order to obtain the asymptotic geometry. It is sufficient to note that for 
cylindrical symmetry the flux tube has the property of Tq — T z z , where TJ? is the energy-momentum 
tensor. This means that the solution will have a symmetry under boosts along the string axis, i.e., a = c. 



The Kasner conditions (1.2) then leave only two options; either the locally flat case: 

a = c = , 6=1 (1.4) 
which we shall refer to as the cosmic string branch, or 

a = c = 2/3 , 6 = -1/3 (1.5) 



which is the same behavior as that of the Melvin solution |19|. We shall therefore refer to (1.5) as the 
Melvin branch. 

The Melvin branch does not provide the same characteristics used to describe a "standard" cosmic 
string, and thus it has been disregarded in previous investigations. However, the Melvin-like solution 
seems completely well-behaved and one may wonder whether it exists at all (asymptotically), or maybe 
it is just an artifact of the too general reasoning above. Although several authors |l7], [l8[ [2C], ^l], ^2) 
explicitly mention this possibility of a second Melvin branch of solutions in the Abelian Higgs system, 
it has never been properly discussed. 

In this paper we show that indeed the Melvin branch actually exists in the Abelian Higgs model. It 
turns out that each conic cosmic string solution has a "shadow" in the form of a corresponding solution 
in the Melvin branch. These two solutions, conic cosmic string and asymptotic Melvin, are found only in 
a part of the two dimensional parameter space of the system, bounded by the curve of maximal angular 
deficit of 2tt. 

There also seem to be some open questions in the literature |l(| [2(| about the nature of the solutions 
beyond the border of maximal angular deficit of 2ir. We will see that also in this region two types of 
solutions coexist, both of them with a finite radial extension. 



2 General analysis of the Abelian Higgs system 

The action of the gravitating Abelian Higgs model is: 

S = J d A xV\J] (^-D^*D^ ^($*$ - v 2 ) 2 -F^ + j^T^j (2.1) 

where 1Z is the Ricci scalar, F^ u the Abelian field strength, $ is a complex scalar field with vacuum 

expectation value v and = V M — ieA^ is the usual gauge covariant derivative. We use units in which 

n = c = i. 

Because of the cylindrical symmetry of the source, we will use a line element of the form: 

ds 2 = N 2 (r)dt 2 - dr 2 - L 2 {r)d(j) 2 - K 2 (r)dz 2 (2.2) 
and the usual Nielsen-Olesen ansatz for the +1 flux unit: 

$ = u/(r)e # , = -(1 - P{r))d(f> (2.3) 
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This gives rise to the following field equations for the Abelian Higgs flux tube: 

(™W + f^(l-f)-^/ = ,2.4) 



NKL \ v L 2 
L ( NK 



NK V L 



P' - e 2 v 2 f 2 P = (2.5) 



With the line element (2.2), the components of the Ricci tensor are: 

(LKNJ v r__N^_£^_Kl 
NLK ' r N L K 

V J- ( NKL 'y n *- ( NLK 'y (2 6) 

K *~ NLK ' NLK [2 - b) 

The source is described by the energy-momentum tensor with the following components: 

To = P = £s+£v + £ w + U 

1~^ r — p r — ■ S s ~\~ Eyj ~\~ 11 

^0 — P(p —~ E s E v ' E w ~\~ U 

V = - Pz =p (2.7) 



where: 

v 2 P' 2 v 2 P 2 f 2 A?; 4 



It turns out to be convenient to use Einstein equations in the form (T = 7^ ) 



7V = -8wG(T^ - ~ 5 „„T) (2.9) 



from which one obtains: 

(LKN')' 
NLK 

(NKL')' 
NLK 

(LNK')' 
NLK 

and instead of the "radial" part of ( |2.9| ), we take the following combination: 

N' L' V K' K' N' , s , 

~N~L + T ~K + ~K ~N = Pr = (£s + £v ~ £w ~ M) (2 ' 13) 

which is not an independent equation but serves as a constraint. In vacuum, the right-hand-sides of 
these equations vanish, and the first three of them are trivially integrated. In this way we may get 



= 4irG(p + p r +v<j>+Pz) =8ttG(s v -u) (2.10) 
-4irG(p -Pr+P<p- p z ) = -8irG(e v + 2e w + u) (2.11) 
= -AttG( P - Pr - p^ + Pz ) = 8ttG(£„-u) (2.12) 



back Kasner's line element (1.1). This is therefore the asymptotic form of the metric tensor around any 
(transversally) localized source and especially around an Abelian Higgs flux tube. Moreover, it is easy 
to get convinced that due to the symmetry under boosts along the string axis, K = N. 

The equations become more transparent if we express all lengths in terms of the scalar characteristic 
length scale 1/VXv 2 (the "correlation length" in the superconductivity terminology). We therefore 
change to the dimensionless length coordinate x = \J Xv 2 r and we introduce the metric component 
L(x) = V \v 2 L(r). We also introduce the two parameters a = e 2 /X and 7 = 8irGv 2 . In terms of these 
new quantities, we get a two parameter system of four coupled non-linear ordinary differential equations 
(the prime now denotes d/dx): 
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(N 2 L/'Y , ( 1 _f2_P"\j = (2 , 4) 



N 2 L V L 2 
L fN 2 P' 



N 2 V L 



-af z P = (2.15) 



(L/V/V')' .Y P' 2 I,, , 2 , 2 



n^-^-^j (2 - i6) 

^ J* + *£ + hi-M (2.17) 



7V 2 L ' V2ai 2 i 2 4 

We have also to keep in mind the existence of the constraint ( |2.13 ), which gets the following form 



N' ( V N'\ ( f 2 P' 2 P 2 f 2 1, 9n9 . 



N \ L N J V 2 2 " L2 2i2 4 

In order to get string-like solutions, the scalar and gauge fields should satisfy the following boundary 
conditions: 

/(0) = , lim f(x) = 1 

P(0) = 1 , lim P(x) = (2.19) 

X — >oc 

Moreover, regularity of the geometry on the symmetry axis x = will be guaranteed by the "initial 
conditions" : 

L(0) = , 7/(0) = 1 

7V(0) = 1 , 7V'(0) = (2.20) 

The purpose of the present paper is to map the two dimensional a-j parameter space, and thereby 
to classify all the string- like solutions of the Abelian Higgs system. It is well-known that, even in fiat 
space, the field equations can only be solved numerically. However, much can be said about the nature 
of the solutions even without explicitly solving the field equations [|[ [l6], |l7], glj . 

The standard Abelian Higgs string solution, usually considered in the literature, has a vanishing 
gravitational mass. This means that the spacetime around the string is locally flat except in the core 
of the string, while there is a non-trivial global effect namely a conical structure of the space which is 
quantified by an angular deficit. However, this does not at all saturate all the possibilities of solutions 



of (2.14 )-(J2.17 ). There are further types of solutions with the same boundary and "initial" conditions, 
(j2.19D-(|2.20D which are not asymptotically flat but have interesting physical interpretations. In this 
paper, we will show that a point in the a-7 plane always represents two solutions, except at the curve 
representing angular deficit of 2tt. The various solutions are distinguished by, among other things, their 
asymptotic geometries. 

For analysing the solutions and obtaining the above-mentioned features and some additional ones, 
we introduce the Tolman mass (per unit length), M: 

roc poo / p/2 1 \ 

GM = 2-kG J drN 2 L(p + p r +p^+p z ) = jJ^ dx N 2 L 1^— - -(1 - / 2 ) 2 J (2.21) 
Using the field equations, one can show that J18[: 



GM = - lim (LNN') (2.22) 
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X — >oc 



We also define the magnetic field B: 

B = 



1 dP(r) 7 / e \ P'(x) 



eL{r) dr a \8ttGJ L(x) 



(is) W < 2 - 23 > 
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and the dimensionless parameter B = 8ttGB(0) / 'e. We find that the central value of the magnetic field 
(its value in the core of the string) can be expressed as p8| : 



B = 1 + 2GM - lim (N 2 L' 



(2.24) 



The asymptotic form of the metric tensor is easily found by direct integration of the two Einstein 



equations (2.16 )-( 2.17) using the boundary conditions and the definitions of M and B. It is of the 
Kasner form 



N(x) = K{x) 



KX 



L(x) ~ (3x b 



with: 



2GM 



6GM + 1-B 



b = 



2GM + 1-B 
6GM + 1-B 



k 2 [3 = 6GM + 1-B 



(2.25) 



(2.26) 



The constant k appears free in the asymptotic solution, but it is uniquely fixed in the complete one by 
the bounda ry co nditions on the metric. 

D ue t o ( [2.18 ) we get the following relation which is equivalent to the quadratic Kasner condition in 
Eq. 0: 

GM(3GM + 1 - B) = (2.27) 

This immidiately points to the possibility of two branches of solutions corresponding to the vanishing 
of either factor in Eq. ( [2.27 ). These two possibilities of course represent the two branches already 
discussed in Section 1. 



Consider first the cosmic string branch (a = c = 0, b = 1) in a little more detail. In this case, we 
find a simple physical interpretation to the parameters in Eq. ( |2.25 ). the constant /3 = L'(oo) defines 
the deficit angle 8<fr = 2ir(l — (3), while the constant n = N(oo) is the red/blue shift of time between 
infinity and the string core. Moreover, relations (2.22), ( [2.24 ) give: 



M = 



B= 1-, 



1 - 



2tt 



(2.28) 



(2.29) 



That is to say, the Tolman mass vanishes and the central magnetic field is directly expressed in terms 
of the red/blue shift k and the deficit angle S<fi. 



Then consider the Melvin branch (a = c = 2/3 , b = —1/3). In this case there is no simple physical 
interpretation of the parame ters (3 and k in Eq. ( 2.25 ). As for the Tolman mass, we notice that it is 
non-zero, but the equations ( 2.22 ), ( |2.24| ) lead to a simple relation to the central magnetic field fl8|| : 



B=l + 3GM 



(2.30) 



For more discussion of these general relations (and some other ones), we refer to Ref. |Ts)| . In 
the remaining sections of this paper, we shall construct explicitly the various types of solutions to 
J2.14 )-(2.17) by using a relaxation procedure to solve the four coupled differential equations numerically. 
The constraint (2.18) has been used for estimating the numerical errors. 



3 Open solutions: Cosmic strings and Melvin branch 

The existence of cosmic string solutions in the Abelian Higgs system is very well established. However, 



the existence of the solutions of the second type, i.e. the Melvin branch type which is implied by (1.5) 



and (2.27), has not been properly studied. We have found that for any conic cosmic string solution, 
an associated solution exists in the Melvin branch. The main difference with respect to the cosmic 
string branch is the fact that asymptotically the azimuthal circles have vanishing circumference. A 
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related difference is the non- vanishing total mass (per unit length) of the Melvin-like solutions. Figure 
|l| shows an example of the conic and the Melvin-like solution at the point (a, 7) = (2, 1.8) in parameter 
space. Only the (square root of the) metric components, N and L, are plotted, as the scalar and vector 
fields for the various solutions deviate only very little from the standard and well examined cosmic 
string configuration. The two branches are obviously quite different. The cosmic string spacetime is 
asymptotically flat as N is constant (actually N = 1 here as a = 2 is the Bogomol'nyi limit |23|, 53) 
and L becomes proportional to x far from the core. In contrast hereto, the associated solution in the 
Melvin branch has N oc x 2 / 3 and L oc x^ 1 / 3 in the asymptotic region. Therefore the circumference of 
circles lying in planes perpendicular to the core with x — as center will eventually decrease as one 
moves to larger x. 

The curves shown in Figure |l| are representative for the cosmic string and Melvin branches. In the 
first case, increasing 7 at fixed a will simply shift N towards a higher constant value in the asymptotic 
region, whereas L will become less steep (meaning larger angular deficit). An increase of a when 7 is 
held constant will have the opposite effect, i.e. shifting L to lower asymptotic values and decreasing 
the angular deficit. This observation holds for all solutions (including the cylindrical solutions and the 
closed solutions to be described in Section 0): only the 7/a ratio seems to matter for the asymptotic 
behavior. Still, because of the region near the core, no obvious symmetry is observed that could 
render the parameter space effectively one dimensional. In the case of the Melvin branch, an increased 
7 (again at fixed a) tends to flatten N, whereas the global maximum of L takes a larger value at higher x. 

Moving towards more massive strings in the parameter space (larger 7 to a ratio) the angular 
deficit of the conic solution will increase until it reaches a maximum of 27r. The corresponding solution 
represents an asymptotically cylindrical manifold with N and L both asymptotically constant. One 
such solution exists for any a and the curve 7* (a), where 8<j> = 2n, plays a very special role in the 
classification of the various solutions. First of all, it obeys a power law jl(|. Secondly, as one approaches 
a point on this curve in parameter space along the Melvin branch the solution converges to the same 
cylindrical solution as for the cosmic string; that is, the conic solutions and the Melvin-like solutions 
coincide in a single asymptotically cylindrical solution on the curve of maximal angular deficit. 

Figure ^ illustrates how the central magnetic field B varies throughout the parameter space for both 
the cosmic string and the Melvin-like universes. Notice how the two surfaces intersect along the curve 
of maximal angular deficit, where the two branches coincide in the cylindrical solutions. 



4 Closed solutions: Inverted cone and singular Kasner solu- 
tions 

The issue of "supermassive cosmic strings" in the Abelian Higgs system was discussed at early days 
fl6[ po[ . Two types of closed solutions were found. One |l6| in which N(x) vanishes at a finite distance 
from the axis (x — x max ) while L(x) diverges at that point. In the other, N(x) stays finite for any x 
while L(x) decreases linearly outside the core of the string, and vanishes at a finite distance from the 
axis. The geometry is still conic but of an inverted one where the apex of the cone is at the point where 
L = Jill ||(| . Using the same numerical methods as in the previous section, we have found that these 
two types of solutions are encountered just as the curve 5<f> = 2tt is traversed. But the lack of asymptotic 



behavior for closed solutions forces us to replace the boundary conditions (2.19) by the condition 



/(0)=0 , f(x max ) = l 

P(0) = 1 , P(x max ) = (4.1) 

This ensures that the solutions are unit flux tubes, and makes the numerical work simpler. For quite 
large radial extension this technicality is not of any relevance as the variables / and P describing the 
scalar and vector fields, respectively, converge quite fast to their values at "infinity". 

Figure || shows an example of both kinds of solutions at the point (0,7) = (2,2.05) in parameter 
space. This is just above the curve of maximal angular deficit, which lies at 7* (2) = 2 for fixed a = 2, 
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and the solutions still have considerable sizes. Strictly speaking, these two types of solutions have 
no asymptotic behavior as they are both representing closed solutions that pinch off at a finite radial 
extension (x max ). But for slightly supermassive configurations we have 

N(x) = K[x) ~ n(x max - x) a , L(x) ~ I3(x max - x) b (4.2) 

The inverted cone behaves according to Eq. (|0|) and the singular Kasner solution behaves according 



toEq. flOp . 

The curves shown in Figure [| are representative for the solutions above the curve of maximal 
angular deficit. As 7 increases (a fixed) both solutions become smaller in radial extent. For the 
inverted cone N gets shifted to a lower constant value and L intersects zero at smaller x. Likewise 
Xmax decreases for the singular Kasner solution causing L to diverge at similarly low x. We have found 
that for any given a and 7 > 7* (a), the solution with geometry as an inverted cone has the larger size. 



As mentioned, the solutions shrink when 7 is increased (a fixed). At some point the spacetime 
described by the solution becomes smaller than a few times the characteristic core thickness which 
scales as a -1 / 2 |(|. At even higher values of 7 the choice of boundary conditions becomes more and 
more important to fundamental questions as flux quantization and topological stability; the solutions 
exist mathematically but are probably unphysical. 

Figure ^ shows three curves in parameter space which have all been fitted to a power law, 
7 = c\a C2 . The lowest one is simply the curve of angular deficit 2ir. This curve separates the 
two open solutions from the two closed ones, and represents itself a linear family of asymptotically 
cylindrical solutions. The middle one represents the curve in parameter space where the singular 
Kasner-like solutions have a radial extent of 10 a -1 / 2 , i.e. ten times the characteristic thickness of 
the core. Similarly, the upper one represents the curve where the inverted cone solutions pinches off 
at Xmax = 10 or 1 ! 2 . Starting on one of these two curves and moving downwards (towards lower 7 for 
fixed a), the two closed solutions will become more and more similar and will eventually coincide in the 
cylindrical solution. Moving upward instead will yield rather small and increasingly unphysical solutions. 



5 Conclusion 

Cylindrically symmetric configurations of the gravitating Abclian Higgs model have been examined 
by application of relaxation techniques to Einstein equations which simplify in this case to a system 
of coupled ordinary differential equations. Everywhere in parameter space two kinds of solutions 
exist, except at the curve of angular deficit 27r, where only one cylindrical solution exists. The two 
open solutions are the standard asymptotically conic cosmic string and a solution with a Melvin like 
asymptotic behavior. The two closed solutions are the inverted cone and a singular Kasner-like solution. 
Supermassive solutions pinch off very fast and are probably unphysical. 

Here only the unit flux tube, n = 1, has been considered, but introduction of an arbitrary n 7^ is 
expected to give similar families of solutions, which can be examined in an analogous way. 

Finally, it is worthwhile noting that the singularities of static field configurations generally seem to 
be lifted when time-dependence is reintroduced |2^| . Therefore the nature of the singularities of the 
supermassive strings may be understood in the framework of time dependent analysis. 
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Figure 1: Typical examples of the numerical solutions at low 7. Here, at the point (a, 7) = (2, 1.8) in 
parameter space, two solutions coexist. The shown metric components N and L are for the cosmic string 
branch a nd the M clvin branch respectively. Both solutions are open and seem well behaved according 
to Eqs. O) and O). 
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Figure 2: The central magnetic field at the string core plotted for both the cosmic string and the Melvin 
branch. As 7 is increased (for fixed a), the solutions in the two branches will converge towards the same 
cylindrical solution causing the two surfaces to intersect at the curve of angular deficit 5(f> — 2ir. 
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Figure 3: Typical examples of the super-massive solutions. Again two different solutions coexist, here 
at (a, 7) = (2,2.05) which is above the curve of maximal angular deficit. Both the inverted cone and 
the Kasner solutions represent closed solutions, and for the Kasner-like solution L has a singularity at 
Xmax where N = 0. 
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Figure 4: Just above the curve of maximal angular deficit, which represents the cylindrical solutions 
(lower curve), two closed solutions coexist. As 7 is increased (fixed a) the radial extent of these 
solutions decreases. The middle and upper curve represents respectively the singular Kasner-like and 
the inverted cone solutions with sizes of ten times the characteristic core thickness. Above these two 
curves the solutions become unphysical. Notice that the three curves can be fitted to power laws; lower: 
ci w 1.64 A c 2 w 0.275, middle: ci w 1.72 A c 2 w 0.285 and upper: ci w 1.87 A c 2 w 0.330. 
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